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We study analytically and by event-driven molecular dynamics simulations the nonergodic and 
aging properties of force-free cooling granular gases with both constant and velocity-dependent 
(viscoelastic) restitution coefficient e for particle pair collisions. We compare the granular gas 
dynamics with an effective single particle stochastic model based on an underdamped Langevin 
equation with time dependent diffusivity. We find that both models share the same behavior of the 
ensemble mean squared displacement (MSD) and the velocity correlations in the small dissipation 
limit. However, we reveal that the time averaged MSD of granular gas particles significantly differs 
from this effective model due to ballistic correlations for systems with constant e. For velocity- 
dependent e these corrections become weaker at longer times. Qualitatively the reported non-ergodic 
behavior is generic for granular gases with any realistic dependence of e on the impact velocity. 

PACS numbers: 81.05.Rm,89.75.-k,62.23.-c 


Granular materials such as stones, sand, different types 
of powders, or their mixtures are ubiquitous in Nature 
and technology, for instance, in the cosmetic, food, and 
building industries lj. Rarefied granular systems, in 
which the distance between particles exceeds their size, 
are called granular gases [2]. On Earth granular gases 
may be realized by placing granular matter into a con¬ 
tainer with vibrating Q or rotating |d| walls, applying 
electrostatic [f| or magnetic [f| forces, etc. Granular 
gases are common in Space, occurring in protoplanetary 
discs, interstellar clouds, and planetary rings (e.g., Sat¬ 
urn’s) @]. In physics granular gases are of fundamental 
importance: they are direct generalizations of idealgases 
when the restriction of ideal collisions is dropped [2(. 

Ergodicity is a fundamental concept of statistical me¬ 
chanics. Starting with Boltzmann, the ergodic hypothe¬ 
sis states that long time averages of physical observables 
are identical to their ensemble averages [§]. In the wake 
of modern microscopic techniques such as single parti¬ 
cle tracking [§!] , in which individual trajectories of single 
molecules or submicron tracers are routinely measured, 
knowledge of the ergodic properties of the system is again 
pressing: while time averages are measured in single par¬ 
ticle assays or massive simulations, generally ensemble 
averages are more accessible theoretically. How measured 
time averages can be interpreted in terms of ensemble ap¬ 
proaches and diffusion models is thus imminent [lo|, HU . 

Here we demonstrate in quantitative detail how ergod¬ 
icity is violated even in simple mechanical systems such 
as force-free granular gases. We analytically derive the 
time and ensemble averaged mean squared displacements 
(MSDs) and show that for both constant and viscoelas¬ 
tic restitution coefficients the time average of the MSD 
is fundamentally different from the corresponding ensem¬ 
ble MSD. Moreover, the amplitude of the time averaged 
MSD is shown to be a decaying function of the length of 



FIG. 1: Collisions in a free granular gas with a below-unity 
restitution coefficient lead to a cooling of the gas. Along with 
the reduced kinetic energy of the gas particles, the diffusion 
coefficient in a free cooling granular gas decreases with time. 

the measured trajectory (aging). Comparison to the ef¬ 
fective single particle underdamped scaled Brownian mo¬ 
tion (SBM) demonstrates that this behavior is due to the 
non-stationarity invoked by the time dependence of the 
granular temperature. Ballistic correlations of the gran¬ 
ular gas are shown to be relevant beyond this effective 
SBM description. Our results for generic granular gases 
are relevant both from a fundamental statistical mechan¬ 
ical point of view and for the practical analysis of time 
series of granular gas particles from observations (in par¬ 
ticular, of granular gases in Space) and simulations. 

Granular gas particles collide inelastically and a frac¬ 
tion of their kinetic energy is transformed into heat stored 
in internal degrees of freedom. The dissipative nature 
of granular gases gives rise to many interesting physi¬ 
cal properties Q. In absence of external forces the gas 
gradually cools down. During the first stage of its evolu¬ 
tion, the granular gas is in the homogeneous cooling state 
characterized by uniform density and absence of macro¬ 
scopic fluxes |2|, realized, e.g., in microgravity environ¬ 
ments [121 . Eventually instabilities occur and vertexes 
develop [2|. Here we focus on spatially uniform systems. 

Collisions are quantified by the restitution coefficient 
e = \v' 12 /vi 2 \, the ratio of the relative speed of two granu- 
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lar particles after and before a collision event [l3| ; e = 1 
denotes perfectly elastic collisions, while e = 0 reflects 
the perfectly inelastic case Hi- For 0 < e < 1 the gran¬ 
ular temperature T(t) = m(v 2 )/2 given by the mean 
kinetic energy of particles with mass m continuously de¬ 
creases according to Haff’s law for granular gases [15(, 

T(t) = T 0 /(l+t/T 0 ) 2 . ( 1 ) 

Here Tq - 1 = | (l — e 2 ) r“ 1 (0) is the inverse characteristic 
time of the granular temperature decay, involving the 
initial value of the inverse mean collision time T“ 1 (t) oc 
\/T(t)/m. Weak dissipation (e ~ 1) thus implies r 0 » 
t c . Due to the temperature decrease the self-diffusion 
coefficient D(t) of the gas is time dependent 'i. fT7i| l~p] . 


D(t) = T{t)T v {t)/m = D 0 /(1 + f/T 0 ) 


( 2 ) 


with T v (t) = 3t c (£)/2, Dq = T 0 T v {Q)/m (Fig.Q]). For £ = 
1 we recover normal diffusion with constant diffusivity. 

Most studies of granular gases assume that e is con¬ 
stant. Different approaches consider the relative collision 
speed dependence of e as £(^ 12 ) — 1 — C\v\!^ + 

[19j . The coefficients C\ and C 2 depend on material prop¬ 
erties and the sizes of gas particles [IcJ. The granular 
temperature of the viscoelastic gas scales as T(t) ~ t -5 / 3 
21. 221 implying D(t) ~ £ -5 / 6 171. Especially in mix¬ 
tures of viscoelastic granular particles a crossover from 
superdiffusion to a giant diffusion occurs [ 2 ^, [24| . 

Simulations. We perform event-driven Molecular Dy¬ 
namics (MD) simulations [25j of a gas of hard-sphere 
granular particles of unit mass and radius, colliding with 
constant (e = 0.8, Fig. [2]) and viscoelastic (e (^ 12 ), Fig. [3]) 
restitution coefficients. The particles move freely be¬ 
tween pairwise collisions, while during the instantaneous 
collisions the particle velocities are updated according to 
certain rules: the post-collision velocities are given by the 
velocities before collision and the restitution coefficient e. 
We simulate N = 1000 particles in a three-dimensional 
cubic box with edge length L = 40 and periodic bound¬ 
ary conditions. The volume density is <f> ~ 0.065. 

We evaluate the gas dynamics in terms of the standard 
ensemble MSD ( R 2 {t )) obtained from averaging over all 
gas particles at time t, as well as the time averaged MSD 

(^) = ^A Jo A ([ R (^ + A )- R (0] 2 )^ (3) 


for a time series R(£) of length t as function of the lag 
time A. Here the an gular brackets denote the average 
(<5 2 (A)) = (1/N) Y/,? <5 2 (A) over all N particle traces. 
For an ergodic system, such as an ideal gas with unit 
restitution coefficient corresponding to normal particle 
diffusion, ensemble and time averaged MSDs are equiva¬ 
lent, (R 2 ( A)) = (P{A)} [ISEI. In contrast, several sys¬ 
tems characterized by anomalous diffusion with power- 
law MSD ( R 2 (t )) ~ t a (a ^ 1) or a corresponding loga¬ 
rithmic growth of the MSD, are nonergodic and display 
the disparity (R 2 ( A)) ^ (S 2 ( A)) fiol. Hit [it] iTi7)j |. 
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FIG. 2: Ensemble ((I? 2 (t))) and time averaged ((<5 2 (A))) 
MSDs versus (lag) time (upper graph) and (<5 2 (A)) versus 
length t of the time series (lower graph), from event-driven 
MD simulations of a granular gas with e = 0.8 (symbols). 
While the ensemble MSD crosses over from ballistic motion 
(7? 2 (£)) ~ t 2 for t -C to to the logarithmic law / R 2 (t )) ~ 
log(t) for t to, the time averaged MSD starts ballistically 
and crosses over to the scaling (<5 2 (A)) ~ A/t given by Eq. Q. 


Fig. [2] shows the results of our simulations of a granu¬ 
lar gas with constant e = 0.8. The ensemble MSD shows 
initial ballistic particle motion, (l? 2 (t)) ~t 2 . Eventually, 
the particles start to collide and gradually lose kinetic en¬ 
ergy. The ensemble MSD of the gas in this regime follows 
the logarithmic law (l? 2 (f)) ~ log(f) (red line in F ig. [H 
top) [2] (see also Supplementary Material (SM) 13 1|). 
The time averaged MSD at short lag times A preserves 
the ballistic law (<5 2 (A)) ~ A 2 . At longer lag times, we 
observe the linear growth (S 2 (A)) ~ A (black symbols 
in Fig. [2l top). In addition to this nonergodic behavior, 
the time averaged MSD decreases with increasing length 
t of the recorded trajectory, (<5 2 (A)) ~ 1/t. This highly 
non-stationary behavior is also referred to as aging, the 
dependence of the system on its time of evolution [30(. It 
implies that the system is becoming progressively slower. 
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FIG. 3: MSDs { R 2 (t )) and (<5 2 (A)) as function of (lag) time 
(top) and (5 2 (A)) versus measurement time t (bottom) from 
MD simulations (symbols) of a granular gas with viscoelastic 
e(ui 2 ). We observe the scaling (i? 2 (t)) ~ t 1 / 6 for t to. 
The time averaged MSD slowly changes between the inidi- 
cated slopes (dashed lines). The continuous change of slope 
of {S 2 ( A)) as function of the length t of time traces from slope 
—5/6 to —1 is seen in the inset of the bottom graph. 


We observe the convergence limA-s.t(<5 2 (A)) — (R 2 (t)). 

Fig. [3] depicts the MD simulations results for a granular 
gas with viscoelastic restitution coefficient £(ui 2 ). In this 
case the ensemble MSD scales as (l? 2 (t)) ~ t 1 / 6 for t 3> 
tq. The time averaged MSD does not seem to follow a 
universal scaling law but appears to transiently change 
from the power-law (5 2 ) ~ A'/ 6 at intermediate lag times 
to ( S 2 ) ~ A at longer A, see the bounds derived in SM 
[3l] ] . As function of the length t of particle traces, we 
observe the crossover from (S 2 ) ~ t ~ 5 / 6 to ( S 2 ) ~ 1/t. 

Granular gas with constant e. Let us explore this be¬ 
havior in more detail. The dynamics of a granular gas 
can be mapped to that of a molecular gas by a rescaling 
of time from t to r as dr = \JT(t)/T(0)dt [Hj. Us¬ 
ing Half’s law 0, it follows that r = r 0 log (1 + t/r 0 ). 
As function of r the granular temperature remains con¬ 


stant, and the velocity-velocity correlation function for 
each spatial component decays exponentially [1, 321, 

(u(ti)u(t 2 )) = T 0 /m x exp(—|r 2 - ti|/t„(0)). (4) 

The MSD then has the exponential approach ( R 2 (t )) = 
6D 0 r — 6-D 0 t,,( 0)(1 — exp[— t/t„( 0)]) to the Brownian law 
(. R 2 {t )) = 6 D 0 t [2, 32|. In real time t we find (i 2 > t\) 

{v{ti)v{ti)) = — (1 + U/'To) /3_1 (1 + t 2 /T 0 )~ l3 ~ 1 (5) 


for the velocity correlator (/3 = to/t„(0)). The MSD is 

(R 2 (t)) = 6Dq [to log (1 +t/ro) 

+r v (O)[(l+t/ro)- 0 -l}}. ( 6 ) 

At short times the particles move ballistically, ( R 2 (t )) ~ 
3D 0 t 2 /t v (0), crossing over to the logarithmic growth 
( R 2 (t )) ~ 6D 0 tq log (t/ro) as seen in Fig. [2l top. 

Based on the velocity autocorrelation © we now ob¬ 
tain the time averaged MSD (see details in SM [ST!) 


S 2 (A)^ ~ 6D 0 r 0 A/t 


(7) 


valid in the range To C A C f, where tq is the char¬ 
acteristic decay time of the temperature law ©. This 
result indeed explains the observed behavior of Fig. [2] 
the time averaged MSD scales linearly with the lag time 
and inverse-proportionally with the length t of the time 
traces. Comparison of Eqs. © and 0 demonstrates the 
nonergodicity and aging properties of the gas particles. 

Viscoelastic granular gas. For a velocity-dependent 
£(ui 2 ) the temperature decays like T(t) ~ T 0 (t/ro) -5 ^ 3 , 
and the time transformation reads r = Grg^t 1 / 6 . The 
MSD in this case exhibits the long time scaling (l? 2 (t)) ~ 


‘iSDuT^t 1 / 6 seen in Fig. [0 top. For the time aver¬ 
aged MSD we analytically obtain the bounds (i5 2 (A)) ~ 
A 7 / & Jt and (5 2 (A)) - A/f 5 / 6 , compare the details in the 
SM 31]. These bounds are given by the dashed lines in 
Fig. [3l top. Concurrent to this change of slopes as func¬ 
tion of the lag time, Fig. [0 bottom, shows the change 
of slope of (d 2 (A)) as function of the trajectory length t 
from the slope —5/6 to —1 at a fixed lag time. 

Scaled Brownian motion. For unit restitution coeffi¬ 
cient individual gas particles at long times perform Brow¬ 
nian motion at a fixed temperature defined by the initial 
velocity distribution. For the dissipative granular gases 
considered herein, the temperature scales like T(t) cs 1/t 2 
and ~ 1/t 5 / 3 , respectively. Single-particle stochastic 
processes with power-law time-varying temperature or, 
equivalently, time dependent diffusivity D(t), are well 
known. Such SBM is described in terms of the over¬ 
damped Langevin equation with diffusivity D(t) ~ t“ _1 
for 0 < a < 2 [33|, |34| . SBM is highly non-stationary and 
known to be nonergodic and aging 11.134135|. 
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To study whether SBM provides an effective single¬ 
particle description of the diffusion in dissipative granular 
gases we extend SBM to the underdamped case, 


dv/dt + v/[T v (t)\ = \/2 D(t)/T v (t) x £(f) (8) 

driven by white Gaussian noise £(f) with correlation func¬ 
tion (£i(ti)£j(t 2 )) = —1 2 ) for the components. 

A granular gas with constant e < 1 corresponds to the 
limiting case a = 0, giving rise to the logarithmic form 
( R 2 (t )) ~ D 0 log(t) and the velocity correlation [36] 




T(0)r o (I+^/tq)^ ~ 
tot„(0)(/3 - 1) (1 +t 2 /T 0 )P 


(9) 


This result for the ultraslow SBM formally coincides with 
the velocity correlation function © for granular gases 
in the limit /? 1, in which the velocity correlation 

time t v ( 0) is much smaller than the characteristic decay 
time To of the granular temperature. This is achieved for 
sufficiently small dissipation in the system (e ;$ 1). 

A more careful analysis shows that there exists a dif¬ 
ference of the SBM model to the full result for the gran¬ 
ular gas due to ballistic correlations [3lJ: The restitution 
coefficient affects only the normal component of the ve¬ 
locity of the colliding particles, while the tangential com¬ 
ponents remain unchanged. Therefore, in a granular gas 
with constant restitution coefficient the trajectories of 
particles become more and more aligned, effecting long- 
termed correlations. In the case of viscoelastic particles 
the restitution coefficient tends to unity when the relative 
velocities of colliding particles decrease with time. There¬ 
fore, at long times, the trajectories become more chaotic, 
and the ballistic correlations disappear. In this viscoelas¬ 
tic case ballistic correlations play a significant role only 
for intermediate lag times, while they become suppressed 
at longer lag times, effecting a gradual crossover between 
the bounding behaviors (|<5 2 (A)^ ~ A 7 / 6 /t and ~ A/t 5 / 6 

(Fig. [3] top). In contrast, for a constant e the ballistic 
correlations are relevant during the entire evolution of 
the gas, canceling the leading term of the SBM result for 
the time averaged MSD [31 1. 

Conclusions. The occurrence of nonergodicity in the 
form of the disparity between (long) time and ensem¬ 
ble averages of physical observables and aging, is not 
surprising in strongly disordered systems described by 
the prominent class of continuous time random walk 
models involving divergent time scales of the dynamics 
Em EH 0 113. Examples include diffusive motion in 
amorphous semiconductors, structured disordered envi¬ 
ronments, or living biological cells 11 j. 

In contrast to such complex systems, we here demon¬ 
strated how nonergodicity arises in force-free granular 
gases. This may seem surprising for such simple mecha¬ 
nistic systems. However, the physical reason for the non¬ 
ergodicity is due to the strong non-stationarity brought 
about by the continuous decay of the gas temperature. 


For a constant restitution coefficient the MSD ( R 2 (t )) 
grows logarithmically, while the time averaged MSD 
(c> 2 (A)) is linear in the lag time and decays inverse pro¬ 
portionally with the length of the analyzed time traces 
(aging). We derived the observed nonergodicity and the 
aging behavior of granular gases from the velocity au¬ 
tocorrelation functions. We note that aging in the ho¬ 
mogeneous cooling state of granular gases was reported 
previously (37], however, it was not put in context with 
the diffusive dynamics of gas particles. 

The decaying temperature of the dissipative force-free 
granular gas corresponds to an increase of the time span 
between successive collisions of gas particles, a feature 
directly built into the SBM model [34[. As we showed 
here, SBM and its ultraslow extension with the logarith¬ 
mic growth of the MSD indeed captures certain features 
of the observed motion and may serve as an effective 
single-particle model for the granular gas, which is par¬ 
ticularly useful when more complex situations are con¬ 
sidered, such as the presence of external force fields. 

Granular gases represent a fundamental physical sys¬ 
tem in statistical mechanics, extending the ideal gas 
model to include dissipation on particle collisions. Gran¬ 
ular gases are a reference model in granular matter 
physics [381 ] with applications ranging from interstellar 
clouds or planetary rings to technologies in food and con¬ 
struction industries [2, l38|. Our results shed new light on 
the physics of granular gases with respect to their viola¬ 
tion of ergodicity in the Boltzmann sense. Moreover, we 
obtained the detailed influence of ballistic correlations in 
the granular gas dynamics. Both results are important 
for a better understanding of dissipation in free gases as 
well as the analysis of experimental observations and MD 
studies of granular gases. 
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In this Supplemental Material (SM) we present the details of the derivation of the results 
presented in the main manuscript as well as an additional Figure comparing the results for 
a viscoelastic granular gas with the numerical evaluation of the time averaged mean squared 
displacement. 


GRANULAR GAS: CONSTANT RESTITUTION COEFFICIENT 


The time-averaged mean-squared displacement (MSD) for the granular gas with constant restitution coefficient, see 
Eq. (3) in the main text, may be presented in the following form 


(<5 2(A) ) = T ±^J* * dt ' [{R 2 {t' + A)) - < R\t')) - 2A (t\ A)] 
where the MSD (i? 2 (t)) is defined according to Eq. (6) of the main text and 

rt nt-\- A 

A(t, A) =3 / dh / dt 2 (v x (t 1 )v x (t 2 )) 

Jo Jt 


(S10) 


3T o r 2 (0) 


1- I 1 + — 

to 


- 1 + 


T 0 + t 


+ 1 


t + A 
to 


(511) 

This term accounts for position correlations at different time moments t and t + A. It vanishes for underdamped 
scaled Brownian motion (SBM) governed by the Langevin equation as well as for ultraslow SBM (Refs. [34,36] of the 
main text). In the present consideration, this term is non-zero because there are long-term ballistic correlations in 
a granular gas with a constant restitution coefficient e. It arises due to the fact that the normal component of the 
relative velocity of colliding particles decreases, while the tangential one remains unchanged in the course of collisions, 
see Ref. [2] of the main text. Introducing Eq. (1S11I) and Eq. (6) from the main text into Eq. (IS10I) . we get the 
time-averaged MSD in the form 

(512) 


;<S 2 (A)} = ( ( 5 2 (A)) + S(A). 

The first part is equal to the time-averaged MSD obtained for the overdamped SBM 

;«aj) = 6Dm 


t — A Jo 
6D 0 t 0 


“V >„ g ^ +1 ' +A 


t' + T 0 


t - A 


[(t + T 0 ) log (t + T 0 ) - (A + To) log (A + T 0 ) — (t — A + T 0 ) log (t — A + t 0 ) + t 0 logT 0 ] .(S13) 


For r 0 < A < 1 


5 2 (A)} cs 


6 A) To A 
t 


The second part has the form 


2(A) = 


6D o t„(0) 


r t -A 


dt' 


1 


log l A 


A 


(S14) 


t' + T 0 


- 1 


< 0 , 


t — A J o 

where /3 = to/t„(0), see the main text. Introducing the new variable y = A/t! we get in the limit To <C A that 

dy 


; (A) ~ -6D 0 t„(0) ^1 - —I (t, A)^ , where I (t, A) = J^ 


v 2 (! + yf 


(S15) 
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This integral can be taken by parts 


Ht, A) = 4^ (l - £)'’+e log (^) (l - jJ 3+1 +HO (ffl , (S16) 

where 

C(0) = -(J3 + 1)J” dy^^^ 7 +i+V’(/3) , (S17) 

7 = 0.5772... is the Euler’s constant, and ^ (z) = is the so-called digamma function. Finally, we get in the 

limit of t A that 

(^(A)} ~6D 0 t 0 C((3)j. (S18) 

This confirms the linear scaling as function of the lag time of the time averaged MSD and thus the nonergodic nature 
of the granular gas. 


GRANULAR GAS: VELOCITY-DEPENDENT RESTITUTION COEFFICIENT 


Similarly, for the viscoelastic granular gas with e = e{y\ 2 ) the time-averaged MSD may be presented as the sum of 
two parts, see Eq. (IS12I) . The first term corresponds to the time-averaged MSD of the SBM process, which can be 
described by the overdamped Langevin equation, namely 


5 0 2 (A) = 


36 Dqt 0 
t - A 


5/6 r t -A 


dt' 


(t 1 + A) 1/6 — t' 1 / 6 


5/6 


216D 0 < 

7(t —A) 


t ve _ A 7 / 6 - (t - A) 


7/6 


The second term, accounting for ballistic correlations, becomes 


;(a) =6D,x (°) [- 
1 t-A L 


exp -- 


6 T, 


5/6 


,(0) L 


(t 1 + A) 1/6 — t' 1 / 6 


- 1 


This integral can be presented as a sum of three parts 


rt -A 


dt' 


exp - 


6 t. 


5/6 


_0 

r„(0) . 


(t' + A) 1/6 -t ,1/6 


- 1 


ok i A 


ok2A 


ot—A 


> fciA 


' k2 A 


We choose the coefficients k \,2 in the following ranges 

1 <C fci < t 0 A 1 / 5 /t®/ 5 (0) and r 0 A 1 / 5 /^ 6/5 (0) « k 2 < t/A. 
This enables us to evaluate the first integral in Eq. (IS21I) as follows 


o k\ A 


dt' 


exp 


6r; 


5/6 r 


t v { 0 ) - 

The third term in Eq. (IS21I) can be evaluated as 


\t' + A) 1/6 - t' 1 / 6 


- 1 


—fci A. 


(S19) 


(S20) 


[...]dt'+ / [...]<&' + / (S21) 


(S22) 


(S23) 


ot-A 


dt' 


Jk 2 A 

rt -A 

' k 2 A 


exp 


6r, 


5/6 


( i ' + A ) 1/6 - t ' 1 / 6 


a - 6t » 


5/6 


( 0 ) 


77 (0) L 


(i' + A) 1/6 -f' 1 / 6 


- 1 


3617 


5/6 r 


7t„(0) L 


-t 7/6 + (t - A) 7/6 + ((fc 2 + 1) 7/6 - fc 7/6 ) A 7 / 6 . (S24) 















































FIG. S.4: Time averaged MSD (<5 2 (A)) divided by lag time A as function of A from MD simulations (symbols) of a granular 
gas with velocity-dependent restitu tion coefficient. The lines connecting the symbols guide the eye. Red line corresponds 
to numerical calculation of (5 2 ( A)) in Eq. (1S12I) for to = 25, r v = 2, Do = 2. These values ensure the closest agreement 
and are consistent with the parameters of the granular gas as used in MD simulations. Dashed line shows the asymptotic 
(<5 2 (A))/A ~ A 1 / 6 behavior according to Eq. 


Under the assumptions (IS22I) the contribution of the terms (IS23I) can be neglected. Finally, assuming that the second 
term in Eq. (IS21I) is small enough compared to Eqs. (IS24I) for the range of parameters k\$ chosen, we get to leading 
order 

(m) * 36fc 2 1/6 ^ o r 0 5/6 ^ . (S25) 

For larger values of the lag time A, in the range r„(0)f 5//6 /t^ 6 CACI that is opposite to the second condition in 
Eq. (1S22I) we have the upper estimate for the correction S(A) to the time averaged MSD of the SBM process, namely 

|5(A)| < 6Dqt v (0) < (<5g(A)) . (S26) 

Then we get in the limit A Ct that 

(m) - (SM) * • (S27) 

In addition to these analytical estimates, we computed numerically the full expression for the time averaged MSD 
(^S 2 ( A)^ = ^<5q(A)^ + H(A). It agrees well with our MD simulation data, compare the curves in Fig. IS.41 where we 

explicitly plot ^<5 2 ^ /A. It shows that in the range of parameters r 0 , t v and D 0 consistent with the MD simulations 

presented in the main text, the transient scaling behavior (|i5 2 (A)^ ~ A 7 / 6 is realized in a limited range of A. Note 

also that in this range the linear SBM scaling for (^S 2 ( A)^ as prescribed by Eq. (IS27I) is no longer valid. The reason 
is that for large values of A, when A —>• t for any length of the trajectory, in Eq. (1S19I) the evolution of the time 
averaged MSD with the lag time A becomes inherently nonlinear. 






























